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Abstract: Cosmic ray anomalies observed by PAMELA and Fermi-LAT experiments may 
be interpreted by heavy (TeV-scale) dark matter annihilation enhanced by Sommerfeld 
effects mediated by a very light (sub-GeV) U{l)x gauge boson, while the recent direct 
searches from CoGeNT and DAMA/LIBRA experiments may indicate a rather light (~ 7 
GeV) dark matter with weak interaction. Motivated by these apparently different scales, we 
consider a gauge mediated next-to-the minimal supersymmetric standard model (NMSSM) 
entended with a light U(l)x sector plus a heavy sector (Hh, H^), which can provide both a 
light (~ 7 GeV) and a heavy (TeV-scale) dark matter without introducing any ad hoc new 
scale. Through the Yukawa coupling between Hh and the messager fields, the U{l)x gauge 
symmetry is broken around the GeV scale radiatively and a large negative m| is generated 
for the NMSSM singlet S. Furthermore, the small kinetic mixing parameter between 
U(l)x and U(l)y is predicted to be 9 ~ 10 -5 — 10 -6 after integrating out the messengers. 
Such a light dark matter, which can have a normal relic density from the late decay of the 
right-handed sneutrino (assumed to be the ordinary next-to-the lightest supersymmetric 
particle and thermally produced in the early Universe), can serve a good candidate to 
explain the recent CoGeNT and DAMA/LIBRA results. 
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1. Introduction 



The recent indirect dark matter detection experiments like PAMELA |ffl| and Fermi-LAT 
Q found cosmic ray anomalies, which can be interpreted by dark matter annihilation 
or decay (although some astrophysical explanations like pulsars are also possible). This 
inspires the construction of a class of models with a light dark U(l)x sector ||, which 
gives a sub-GeV dark gauge boson. Such a sub-GeV gauge boson plays a key role in the 
dark matter explanation of the cosmic ray anomalies: for the annihilating dark matter it 
can induce large Sommerfeld enhancement and kinetically forbid the hadronic products 
from the annihilation, while for the decaying dark matter it can suppress the hadronic 
activity ||. At the same time, some dark matter direct detection experiments such as 
DAM A/LIBRA @, CDMS II § and CoGeNT § also reported some plausible evidence of 
dark matter, which, including the null result from XENON10 @, may be accommodated 
by a quite light dark matter at GeV scale (~ 7 GeV) with a dark matter nucleon scattering 
cross section a p ~ 10~ 40 cm 2 [O. This has inspired some recent studies on the light dark 



matter JTT|, |TJ]. 

With Sommerfeld enhancement, it seems to us that the dark matter explanation for 
all these experiments must involve three very different scales: the TeV-scale heavy dark 
matter (HDM), the GeV-scale light dark matter (LDM), and the sub-GeV U x (l) dark 
sector. It is then quite challenging to embrace all these aspects in one framework. Firstly, 
it is not a trivial problem to accommodate such a light U(l)x gauge boson at low energy 
without introducing a new scale by hand. As is well known, supersymmetry (SUSY) helps 
to stabilize a scale and, moreover, its breaking usually generates a new scale which is 
encoded in the soft SUSY breaking terms. Thus the crucial task is to obtain a proper 
GeV-scale soft Lagrangian for the Higgs fields in the U{l)x sector. As proposed in S 



and then followed in [13, 14, |15|, [lq| , SUSY breaking (maybe exhibited as soft terms) in 
some hidden sector may be gauge mediated to the U(l)x sector to generate the GeV- 
scale. Secondly, although it is not difficult to construct a GeV-scale U(l)x sector while 
allows for a sub-GeV gauge boson through introducing a very weakly charged Higgs field 
(say Qh9x ~ 0.03), the U(l)x sector with such a light gauge boson will usually also 
predict some other Higgs bosons as light as the gauge boson and the LDM annihilates to 
these bosons very effectively, leading to a very small relic density after freezing out (say 
^■LDMh 2 ~ 10 -4 ). Some studies [12] showed that even with such a small relic density the 
LDM may still generate scattering signals at the dark matter detectors if the LDM-quark 
coupling strength is enhanced enough. Nevertheless, it would be more natural if the LDM 
density is at a normal level (~ 0.1). Thus, the LDM may be understood to be mainly 
produced from the late decay of the ordinary next-to-the lightest sparticle (ONLSP) in the 
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visible sector. This may be a reasonable conjecture since in the presence of some new light 
R— odd state in the U(l)x sector, the ONLSP may decay to this sector with proper time 
scale. 

In this work we try to extend the gauge mediated next-to-the minimal supersymmetric 
standard model (NMSSM) with a light U(l)x sector plus a heavy sector (H^jH^), which 
can provide both a light (~ 7 GeV) and a heavy (TeV-scale) dark matter without introduc- 
ing any ad hoc new scale. In our framework the crucial dynamics is that the HDM couples 
directly through Yukawa couplings with the messenger fields which carry the U (l)x charge. 
The U(l)x gauge symmetry can be broken around the GeV scale radiatively, and a large 
negative m 2 s is generated for the NMSSM singlet S. Interestingly, the small kinetic mixing 
parameter between U(l)x and U(l)y is predicted to be 8 ~ 10 -5 — 10~ 6 after integrating 
out the messenger fields. Such a light dark matter, which can have a normal relic density 
from the late decay of the right-handed sneutrino, can be a good candidate to explain the 
recent CoGeNT and DAMA/LIBRA data. 

This work is organized as follows. In Section II we present the model. In Section III we 
discuss its concrete realization. Finally, discussions and conclusion are given in Section IV. 
In Appendix A, we explain the kinetic mixing and dark-visible interaction. In Appendix 
B, we present the soft terms from HDM-messenger direct couplings. And in Appendix C, 
we give the one- loop renormalization group equations (RGEs) of some soft terms. 



2. Model Building 

Our model is based on the NMSSM with gauge mediated SUSY breaking (GMSB). And it 
has two features: (1) The NMSSM singlet S naturally provides a TeV scale to explain the 
origin of the HDM mass scale; (2) Through radiative correction with 1/167T 2 suppression, 
the GMSB provides a simple way to generate the GeV-scale for the U(l)x dark sector. 
Some previous studies on this line have been carried out [||, 13, 14, [l5|, |16[| . In our study we 



will intensively examine the dark matter phenomenology in the NMSSM extended with the 
U(l)x sector and the extra TeV-scale degree of freedoms, paying special attention to the 
mechanism of the U(l)x breaking at GeV-scale. We find that if the conventional hidden 
sector messengers are slightly charged under U(l)x , then the soft terms in the U(l)x dark 
sector can be at a proper scale. Our work will address the following problems in a coherent 
framework: 

2.1 Generating a Large Negative Soft Mass-Square for S 



As a simple extension of the MSSM, the NMSSM [17| can solve the /i, problem and the little 



hierarchy problem [18], and thus has recently attracted much attention [|H|. However, in 
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the mechanism of the GMSB it is difficult to construct a phenomenologically acceptable 
NMSSM [20, pi | . The key difficulty is that the singlet S couples only to the Higgs doublets 
and thus the soft term m 2 s can not be generated negative enough at the weak scale through 
RGE. To solve this problem, some efforts have been made, e.g., coupling S to extra SU(3)c- 



charged particles [^(J or directly to messengers [22, 23, p4| ]. In our framework, since we 
have extra states (Hh,Hh) which couple to S, we can obtain large enough m 2 s by only 
coupling Hh or directly to messengers (S does not couple to messengers). In fact, this 
is a natural choice since this coupling leads to a large (~TeV 2 ) splitting between the soft 
mass-square nig and fn\{ h at the messenger boundary. That significantly impacts on the 
evolution of the soft mass-square of the dark Higgs field, leading to a negative mass-square 
and breaking the U(l)x in the dark sector. The dynamics of this part is described by the 
superpotential 

Wi = (\SH u H d + |s 3 ) + X h SH h H h 

+ H h (AtTiT 2 + X D D X D 2 ) + X (6,tTiT! + £i, D L>iZ>i + (1 -> 2)) , (2.1) 

where X is the spurion Goldstino field parameterized as X = M+9 2 F, and (Tj, D{) = fi and 
(Ti,Di) = fi form (5,5) representation of SU(5)— GUT group. The matter and messenger 
fields have the assignments under the Z%— symmetry of the NMSSM and the U(l)x ■ 

S -> e™'*S, H h -> e^ 3 H h , H h -> H h , 

[fx] = -[fi] = Qh, [H h ] = -[Hh] = ~Qh h , (2-2) 

while all other fields are neutral under the above symmetries, thus Qf x = Qn h - 
Let us comments on the superpotential: 

(1) The superpotential has a Z% symmetry to keep the HDM stable (the messengers 
(/i, /i) are — odd). According to a recent study [25], the explanation of PAMELA 
through such HDM annihilation with Sommerfeld enhancement is difficult. In par- 
ticular, the maximal Sommerfeld enhanced factor is about 100 for a TeV scale heavy 
dark matter. To explain the PAMELA and Fermi-LAT experiments, for simplicity, 
we assume that the dark matter density in the sub-halo is about three or four times 
larger than the usual value. By the way, to explain PAMELA, we had better re- 
sort to decaying HDM. To let our HDM to decay to dark gauge bosons, we need to 
break the Z 2 symmetry by introducing some new mechanism || . We will not further 
discuss the HDM phenomenology in this work. Instead, we will focus on the LDM 
phenomenology. 

(2) It is important to arrange U(l)x charge to forbid coupling like XfHhfifi, which leads 
to a one-loop tadpole for S in the superpotential after integrating out the messengers: 
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f d 2 9£S with £ ~ Aj-F/167r 2 , which tends to destabilize the weak scale. But at the 
messenger boundary, a large negative mass-square for S is generated at two loop as 

F 2 



m s = ~ ( 3A t + 2A^) Xlj^, (2.3) 

which can be as large as several-hundred GeV, depending on the couplings. For 
example, for M ~ 10 8 GeV at the messenger boundary and taking Yukawa couplings 
as A/j ~ 1, Xt — Ad ~ 0.2, we have m| ~ —(280 GeV) 2 . In this way, it is possible to 
make the NMSSM in the GMSB to have successful electroweak symmetry breaking. 

(3) The Yukawa coupling A^ plays an important role. In additional to generate a large 
m|, a large Xh is also required for having a HDM. From Eq. ( |C.9| ), the one-loop evolu- 
tion of Xh below the messenger scale is approximated as (drop the small contribution 
from A, K and Qh9x) 

^ M »-»Kw^ iog W 1/3 ' {2A) 

We need Xh(M susy ) ~ 0.5 — 1 (depending on the value of v s ) to have a heavy HDM. 
Besides, it makes the HDM to annihilate to some states in the NMSSM effectively 
so that to have small relic density. In this way the HDM can avoid direct detection 
and explain the cosmic ray anomaly by a proper shorter lifetime than the decaying 
HDM with the assumption Q^DMh 2 — 0.12. 

2.2 Generating A Small Negative Soft Mass-Square for the Dark Higgs H 

In our model we assume that the dark sector respects a global SU(N) flavor symmetry 



(it can be gauged to form a non-Abelian dark sector [26, 14], but in this work we do not 
discuss this case). This symmetry is useful because it can protect the light dark matter 
candidate to be stable and allow to construct a simple dark sector without anomaly if we 
require the dark sector has no U(l)x singlet (we will explain why we do not prefer a singlet 
later). The minimal field content includes: (Hi, Hi) carrying U(l)x charge (Qg ,Qhi) and 
forming the (N, N) representation of SU(N); the dark Higgs H carrying U(T)x charge Qh h 
and being a flavor singlet. Then, under the symmetry U(l)x xSU(N), the most general 
superpotential takes a very simple form: 

W dark = XiHHiHi. (2.5) 

The U(l) x anomaly cancellation and the U(l)x neutral condition lead to two equations: 

Q Hh + N(Q Hl +Q Sl ) = 0, Q% h + N(Q H[ +Q%) = 0. (2.6) 
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Note that other U(l)x -charged states in our model are vector-like, and thus do not con- 
tribute to anomaly. Especially, it has a nontrivial solution 



Qh > = ~ (i-^^W^-^) ' (2 ' 7) 



Qh> = [l + ^QH^^^-j ■ (2.8) 

Another solution is trivial, obtained by exchanging the role of Hi and Hi. For any allowed 
N, Qh 1 and Qg t take opposite sign with Qh h ■ This is a required property to assure that 
only H gets negative soft mass-square. 

^From the requirement of a negative m 2 H at the dark scale Qn h is determined to 
take the same sign with Qn h - For pure GMSB, at the messenger boundary, due to anomaly 
cancellation, there is a sum rule for the soft terms: Sx = Tr(Qjm 2 ) = 0, with the trace 
running over all U(l)x -charged fields. But the U(l)x -charged HDM directly couples to 
the messengers in the hidden sector and acquires a large boundary value through Yukawa 
mediation (see Eqs. ( |B.6| ) and ( |B.7| )), which violates this sum rule. Consequently, the non- 
vanishing Sx drastically changes the renormalization of the dark Higgs soft mass-squares, 
driving some of them negative at fin- The trace term is then given by 

Sx = Tr(Q;m 2 ) = Q Hh ( m H h ~ m R h ) 

~ Q»h 7^1)2 l^tiiM) - 5A 2 (M)] ^. (2.9) 

The above estimation is based on the requirement that at the scale M, \t,d 93, A^, 
where g%(M) — 0.9. There is a substantial cancellation between the terms in the bracket 
and thus in the following estimation we set 16g 2 (M) — 5A^ = Ct<7§(-M) with Ct ~ 1. 
Consequently, it generates the low energy value of m 2 H (see Eq. (|C.3| )): 

* ™ 2 „(M) + ^^S X x log g, (2.10) 

where the first term is the small soft term contributed by pure U(l)x gauge mediation. 
Then, we can parameterize the low energy dark Higgs parameter as 

2 / r\ \ 2 



0.16 (9m.\ (9m.) 

V 0.01 / V 0.02 / 
( Qh9x\ ( QH h 9x\ ( At 



-0.2 



V 0.02 J \ 0.01 J V0.2 



GeV 2 . (2.11) 



Here we set F/M = 10 5 GeV, Ct = 0.5 and a typical dark scale [id ~ 10 GeV. With 
a moderate arrangement for Qf 1 gx and At, we readily get > m 2 H (nd) > —1 GeV 2 . 
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Note that we do not need cancellation between the two contributions. In practice, we only 
require that the second term dominates over the first term and takes small value. The soft 
mass-square for (Hi, Hi) can be obtained similarly, which is enhanced by the trace term 
because their U(l)x charge is opposite to H. 

Let us comment on the above charge assignments. First, Qf x gx ~ 0.01 not only 
determines the soft mass scale from U(l)x mediation, but also directly relates with the 
value of 9 discussed later. Next, a small value Qn h 9x ~ 0.01 ont only helps to make the 
HDM to avoid direct detection, but also avoid the unnecessary enhancement by dark gauge 
boson, which will subject to the gamma ray constraint . As for the small Qh9x ~ 0.02, 
controlling the quartic term from D-term, is necessary to generate a larger VEV of the 
dark Higgs, providing the several-GeVs scale for the light dark matter. 

2.3 Predicting a Small Kinetic Mixing Parameter 9 

Since the messengers also carry U(l)x charges, our framework naturally predicts a value 
for the kinetic mixing parameter 9 between U(l)x and Uy(X)'- 

a sr^ 9i9xQyQ x M GUT 

°^ ni 16^ lQg ^~' (2 - 12) 

where nj is the number of fields I that carry hypercharge Q Y and U(l)x charge Q x . At 
first, the contribution from a complete representation of SU(5) cancels exactly due to the 
traceless generators of SU(5). For example, for we have 

2x^xQ /l + 2x f-ij xQ fl )= 2Q /l Tr(T 24 ) = 0, (2.13) 

with T24 a generator of SU (5) that defines hypercharge. In general, the small doublet- 
triplet splittings between messenger fields can be obtained after the SU(5) gauge symmetry 
breaking via the dimension-5 operators Xf^fi/A (with operator coefficient set to be 1) 
where <E> is the 24 representation Higgs field and A is the fundamental scale of the UV theory 
such as string scale M str ing — 3.0 x 10 17 GeV, or reduced Planck scale Mp\ ~ 2.4 x 10 18 



GeV f7|. Note that the SU(5) unification scale is about 2.4 x 10 16 GeV, thus if we take 



A = M s t r i n g, we obtain — £i,t|/£i,.d ~ 0.1. In addition, the 9 parameter is given by 

^2Q /l ^log|^~O(10- 5 ) . (2.14) 

Thus, we require Qf x gx ~ 0.05, which is consistent with the previous parameterization. 
By the way, the RGE effects may also induce the doublet-triplet splitting (see Eqs. (|C.12 ) 



In principle, we can increase QhQx ~ 0.5 and meanwhile set a smaller value At ~ 0.03. It makes the 
HDM still be a candidate for annihilating dark matter with Sommerfeld enhancement. 
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and ( C.13|) ), which could be very small if the corresponding Yukawa couplings are small. 
On the other hand, for large Yukawa couplings, proper splitting is induced even without 
turning back to the high dimension operators. 

2.4 Light Dark Matter Candidate 

Recently, dark matter direct detection experiments showed some hints on light dark matter 



~ 7 GeV. It is natural to relate it with the light U(l)x dark sector [12, [2^]. Although there 
is a small gap between the LDM and dark gauge boson mass scale, it can be explained 
by a small gauge coupling of the dark Higgs (Qh9x ~ 0.02), provided that the Yukawa 
coupling Xi is about 0.5. We will elaborate this problem in the next Section. 

In summary, we depict our dynamics structure in Fig. |]. The hidden sector plays 
a crucial role in our framework: it not only generates all the necessary low energy mass 
scales, but also explains a small 9 in the dark matter phenomenology. 



)M gauge 




Figure 1: The schematic diagram showing our dynamics. Solid and dashed lines denote other 
possible interactions beyond gauge interactions. 



3. Light Dark Matter Phenomenology 
3.1 Vacuum and Spectra of the Theory 

First we check the vacua of the visible sector and the dark sector. For the former, the 
desired electroweak vaccum takes a form v s = (S), v u = (H®) and v<i = {H®). In order to 
ensure the stability of the HDM and avoid the breaking of U(l)x at TeV scale, we must 
have (Hh) = {H^) = 0. We need to be cautious about this because m 2 H (M) and m^> (M) 

h 

are negative and take roughly a value as m|(M) (see Eqs. QB.6| ) and ( JB.7D ). But a negative 
mass-square does not always mean a non-zero VEV in the multi-Higgs system. We can 
prove it by assuming a vacuum with vanishing (Hh) and {Hh), and then check whether such 
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a vacuum leads to a tachyonic direction. In practice, in the complex scalar mass system of 
(Hh, Hfy, the mass matrix reads 



{Xhv s y + m 2 B (X h v s )(Kv s ) + [\ h v s )A\ 



Mjj h =\ v " " • " »h v " °" °' v " " " h . (3.1) 
\(X h v s )(Kv s ) + X h v s A Xh (X h v s ) 2 + m 2 Hh J 

Obviously, it is definitively positive provided that XhV s ~ 2 TeV is much larger than other 
scales in the matrix. This condition can be satisfied for the following reasons. First, 
from the previous parameter estimation, all the soft masses typically lie much below TeV. 
Furthermore, to generate a large v s we require n ~ 0.1 (<S A^). Concretely, the lightest 
boson has a mass-square approximated by 

Ml « (X h v s ) 2 + m\ - ( 2 M ^ fe) ' 2 2 . (3.2) 

n h H h 

where (M^)i2 is the 12-element of Mfj } . This approximation is valid when rn 2 Hh —TUg > 
\(Mjj ) 12 1 - So, in general this is a stable dark matter in the TeV region required by 
PAMELA and Fermi-LAT. 

We further briefly comment on the pattern of the parameter space and the symmetry 
breaking in the NMSSM in our scenario. First, A\ and A K are induced by RGEs, suppressed 
by loop factor. Note that a new contribution A\ h ~ —100 GeV (see Eq. flB,5|) ) affects the 
running of A\ and A K significantly: 

dA\ 

16vr 2 ^ « {2X\A K + 6h 2 t A t + 6g 2 M 2 + ...) , (3.3) 

16 7 r 2 ^^(2X 2 A Xh + 12X 2 A x + ...). (3.4) 

Since is large in our framework, this new contribution is quite sizable. Especially, A K 

(3X 2 +2 A 2 )A 2 rri 

will get a new contribution at order ~ —2-^ — jy^zjr h jj l°g m ~ ^(10) GeV. As a result, 
in general one can not expect a very light R— axion (CP-odd) a in the spectrum. But in 
case of small n and A, and v s ^> v, some parameter space still allows for m a < Im^ and 



consequently the R— axion solution to the fine-tuning problem may be accommodated ]18 

4 



Anyway, these trilinear terms are small compared with m%, and thus the electroweak and 



Z$ breaking in the NMSSM is dominantly driven by the negative m|. Approximately, we 
get a ti s » t) limit through 

Vs ~ OH „ 0(5)TeV, (3.5) 

K 

which is readily achieved by a small k ~ 0.1. Incidentally, a small k is a safe choice to 
stabilize the HDM mass scale and moreover is favored by keeping the theory perturbative 
up to the GUT scale. In conclusion, with the effects of the Yukawa coupling between the 
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HDM and the messengers, our scenario is capable of providing a proper solution to the 
NMSSM. 

The U(l)x symmetry breaking and the spectrum in the dark sector can be analytically 
studied. The total scalar potential is V = Vr> + Vf + V so ft, with each term given by 2 

V D =\g 2 x [{Q Hl \H l \ 2 + QH l \H l \ 2 + QH h \H\ 2 )+ix]\ (3.6) 
V F = \\imml 2 + |A,| 2 |tf | 2 (\Hi\ 2 + \H\ 2 ) , (3.7) 
V S oft =m 2 Hl \Hi\ 2 + m 2 Ri \H l \ 2 + m 2 H \H\ 2 + {\ l A h HH l H l + h.c) . (3.8) 

Among the soft terms, m 2 ^ and fn 2 Hi are positive while m 2 H is negative. A\ t is purely RGE 
induced, roughly given by (see Eq. (|C.6|) ) 

A ^ - ~ TTT^iQmgxfmxlog — ~ -10" 2 GeV, (3.9) 

which is much smaller than the typical scale in the dark sector and thus is not a relevant 
soft parameter although it controls the mixing between Hi and Hf, H is the Higgs field 
which breaks U(l)x gauge symmetry. Its potential is simply a complex 4 theory, where 
the negative m 2 H and quartic term from D— term stabilizes the potential at the minimum 

(H) = v H = ~ O(10) GeV. (3.10) 

Qh9x 

The dark spectrum can be at the required several-GeV scale simply by setting A; ~ 0.5. The 
dark gauge boson mass is given by rax = V%Qh9xvh = V^\mn\ — 0.2 GeV, depending 
only on the negative Higgs parameter. To calculate the dark spectrum and the interactions 
in the dark sector, we take a unitary gauge to eliminate the Goldstone boson from the 
spectrum and write 

H = v H + ^. (3.11) 

The CP-even state Hr does not mix with other states and gets its mass from the quartic 
term {D— term). Since the D— term is determined by gauge coupling, at tree- level Hr is 
exactly as light as the dark gauge boson. The LDM can annihilate into such light bosons 
too effectively and thus the freeze-out relic density is too low, which will be discussed in 
the next section. 

Now we study the states from the superfields (Hi ,H[). The mass-square matrix of the 
complex scalars in the basis of (H^,H{) is given by 

Xfv 2 H + m 2 Hi XiA Xl v H 
XiA Xl v H X 2 v 2 H + m 2 Si 

2 The study in |Ie| observes that the effective Fl-term £i oc 8 generated by the mixing D— term between 
U{l)y and U(l)x is able to generate the proper U(l)x breaking, which is ignored in our study because 
is small. 



Mf = \ 'r A ^ Q ). (3.i2) 
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Here we do not consider CP-violation and thus all parameters are taken as real. The soft 
trilinear term is only a RGE effect and is a small perturbation to the diagonal elements. So 
the two mass eigenstates are approximately same as the interaction states (H[,Hi), with 
the eigenvalues given by the two diagonal elements in M 2 . Because of the positivity of 
Tn 2 Hi and m 2 ^ , these bosons are heavier than the Dirac fermion formed by the two Weyl 
fermions % = (Hi, Hi), whose mass is given by Ml = XiVfj- Such a U(l)x -charged Dirac 
fermion \ serves as the LDM candidate. 

Note that in our framework we do not choose a singlet type dark sector as in [15], where 



the dark sector superpotential is NH'H, with N being a singlet and Qh> = —Qn h , and 
no SU(N) flavor symmetry is introduced. Since in that case N is a singlet with m 2 N < 
at obtained from renormalization only, the LDM is always a singlet-like scalar and its 
couplings with quarks are suppressed by an extra mixing factor 1/5^, with 1/5a measuring 
the fraction of the charged component H 1 in the LDM, given by 

\MA Xi vh\ hgx 

In this estimation we assumed an ideal case that the pure U(l)x mediation contribution 
to mjj and m 2 ^ equals to the renormalization contribution from Sx- So in that scenario, 
besides the suppression from 9 2 ~ 10~ 10 , the LDM-nucleon scattering will be further 
suppressed by a factor 1/S\ ~ 10 -8 , rendering the cross section unacceptably small. 

Finally, we comment on H (fermionic component of H) and dark gaugino X. They 
have a Dirac mass mux = m x — 0.2 GeV. And X also has a heavier Major ana mass term 

~ 0.5 GeV for the choice Q^gx — 0.01. This will lead to a seesaw-like spectrum, i.e., 
the lighter one is very light, even as light as tens of MeV. Provided that the SUSY breaking 
scale is high enough, saying yF ~ 10 9 GeV, this particle will be the LSP (otherwise, we 
have to make sure that after decoupling, it decays away before the BBN). 

3.2 Light Dark Matter Relic Density 

In this paper we mainly discuss the LDM phenomenology and try to explain the CoGeNT 
and DAMA/LIBRA results together with other null results from XENON100, XENON10 
and CDMS(Si). We use the latest data analysis in [{Tcfl , which showed that the combination 
of DAMA/LIBRA and CoGeNT data can be well accommodated by a LDM with a mass 
of ~ 7 GeV and an elastic scattering cross section with the nucleon of ~ 2 x 10 -40 cm 2 . 
Moreover, it showed that such a LDM is not excluded by other null results. 

In our framework we have such a LDM from the dark sector, the Dirac fermion \. 
However, a proper relic density for this LDM is hard to obtain from the standard freeze- 
out thermal production. In fact, there are two annihilation channels for this LDM: one is 
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directly to the dark gauge boson with a rate oc 47r(QH i 5x) 4 /"i| and the other is to Hr 
with a rate oc Ai:Xj/m 2 L . Clearly, without the suppression of any large mass scale (e.g., 
a weak scale heavy field in the propagator), the only way to keep the LDM as a weakly 
interacting massive particle (WIMP) with a typical weak reaction rate ctq ~ 10~ 8 GeV -2 
is to set Qn t <?x ; A/ i$ 0.03. But such a smaller A/ implies vh must take several hundred 
GeVs to keep the mass of x ls about 7 GeV. Anyway, in principle it is a viable solution, 
for example, by taking 

Qh9x ^ 0.001, A; ~ 0.03, 9 ~ 10" 4 , (3.14) 



and keep mx as light as 0.2 GeV, i.e., m 2 H (nd) — —0.04 GeV 2 . But from Eq. ( p. lip we 



have to choose Qf x gx — Qh h 9x — 0.2. Although this solution has a virtue that it allows 
the HDM to be the Sommerfeld type instead of decaying HDM (since HDM couples to dark 
gauge boson with a large strength), it is at the price of a surprisingly large U{l)x charge 
hierarchy between the different fields, e.g., Qn x : Qh h = 1 : 200. So we propose that 
the LDM abundance is produced by the late-decay of the ordinary next-to-the lightest 
supersymmetric particle (ONLSP) in the NMSSM (the collider phenomenology of such 
ONLSP decay to dark states is studied in [14]). In this solution, A/ takes a large value so 



that the LDM annihilates to Hr very fast and eventually leaves a small abundance after 
decoupling. However, in the presence of a light U(l)x sector, the ONLSP will dominantly 
decay to the U(l)x -charged dark states. If this decay happens after the decoupling of the 
LDM (typically ~ 10~ 5 s), then the number density of the ONLSP is transferred to the 
LDM. Of course, a rather large number density is needed because the relic energy density 
of the LDM is proportional to its mass. But if the ONLSP has a very weak annihilation, its 
relic number density can be quite large. In the following we discuss this issue quantitatively. 

First we consider the lightest neutralino N\ as the ONLSP. If N% dominantly decays 
to dark Higgsino and dark Higgs, its lifetime is estimated to be 

where /g is the fraction of bino in N±. In order for this decay to be late enough, the bino 
component should be highly suppressed ~ 10 -4 . 

Then we assume the right-handed snuetrino (sRHN) as the ONLSP. Such a sRHN is 
present in the NMSSM extended with a right-handed neutrino, which was used to explain 
the light neutrino masses by seesaw mechanism p9| , |30[ . We consider a simple model 
with only one flavor RHN (denoted as N) and lepton doublet introduced. Its relevant 
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superpotential is given by 
W N = 



Y N LH U N + ^N 2 + ^SN 2 + »H u H d , 



(3.16) 



2 2 

where \x = Xv s . Depending on the ^-charge assignment, Mn or \$n can be turned off. In 
the following we focus on the case with A$n = 0. Further, in GMSB the soft terms involving 
the SM singlet N can be dropped safely because they are all generated by RGE effects from 
the coupling to L and H u (such renormalization effect is suppressed by Y N ~ 10~ 5 in the 
low-scale seesaw). The LR- mixing is naturally suppressed. In the CP-eigenstate basis of 
sleptons (V*y. , Nl , v*_ , ) , the mass matrix is 

(m 2 ~+D 2 F 2 + m D M N \ 


rril + D 2 F 2 -m D M N 

B N M N 



Ml 



m % + M N + b nM n 



\ 



m 2 ~+M 2 N 



(3.17) 



/ 



where D 2 = 0.5m| cos(2/3) and the mixing parameters F 2 



-/iitid cot /3 



state v\ with mass-square rriy is dominated by N*, provided that 



m% + M 2 - B N M N ^M 2 <m 2 ~ + D 2 . 



The lightest 



(3.18) 



We have used the fact that the splitting is small, so the mass eigenvalues are nearly the 
four diagonal elements. Concretely, the component of v*_ is given by v\ D C^[v*_ with 

F 2 - m D M N F 2 - m D M N 



CI 



5m 2 2 



(3.19) 



where 5m 2 2 * s the mass-square splitting between the two mass eigenstates of 3-4 block in 
the matrix Eq. ( 3.17] ). Depending on the mass splitting and M^v, the fraction covers over 
a wide region: 



m v M% 
5m\ 2 



10" 



10" 



(3.20) 



where we used the seesaw formula for the light neutrino mass scale m u = m 2 D /Mi^ ~ 0.1 
eV. 

As the ONLSP, the D\ has two decay channels to dark sector through its left-handed 
slepton component. One is the interesting three-body decay via v\ — > vl+Hi+Hi mediated 
by bino, as shown in Fig. [2|, and the decay lifetime is [14 

-i 



~2.6 x 10" d s x 



16vr 2 P(m ?1 /M : 



10" 



10" 



10" 



300 GeV 



9 J Q 2 a x m Vl P(m v jM 1 



(3.21) 
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The other channel is v\ — > vl + X (also see Fig. ^), but is suppressed by an additional 
helicity factor (m^/Mi) 2 and typically several times smaller than the three-body decay 



1 14]. Moreover, it can also the decay into Goldstino i>\ — >• ul + G, as shown in Fig. |2j, 
which is suppressed by the SUSY-breaking scale y/F > 10 3 TeV. So, the v\ ONLSP mainly 



decays to the U(l)x charged dark states before BBN era (> Is) and thus can provide a 
proper LDM density. 



$19x9 
Hi 

Figure 2: The v\ decays to dark states and gravitino 






To end up this section, we point out one merit of the LDM from late decay. In Ref. [32] 
it was shown that if the LDM reaches its relic density via annihilating to SM fermions, 
then the required LDM-nucleon scattering cross section generally implies antiproton excess, 
leading to some tension. But obviously our LDM scenario evades this constraint. 

3.3 Explanation of CoGeNT and DAMA /LIBRA Results 

The study in [|l0| suggests a ~ 7 GeV light dark matter with an elastic scattering cross 
section with the nucleon of ~ 2 x 10 -40 cm 2 . In the following we study the LDM interaction 
with the nucleon. 

Using the method described in f33|| , we derive the effective interaction between the 
LDM and the nucleus. The microscopic interaction is presented in Appendix. [A|. Due to 
the kinetic mixing, the LDM interacts with quarks, mediated by the dark gauge boson. It 
is the basis of the effective theory describing the LDM-nucleon interaction. This effective 
theory is obtained by calculating the quark and gluon operators in a nucleon state, such 
as (nlfj^fln). Then we obtain an interaction: 

C n J c = J£ (X^ + QswZp) + fnX^n + ffX&fp, 

J x = + Qsjgxx-fx + \(QH t - Qa^xx-f-ysx 

+ iQigx [d^HtHj - d^HjH^j + • • • , (3.22) 
where the dots denote the irrelevant contributions like the gauginos. From Eq. ([^7]) we get 



Qh { +Qhi = Qn h ■ And the term proportional to \Q Hl -Qfii I = \J ^m~\QH h \ < \Qh 1 +Qh 1 I 
induces spin-dependent scattering. But it is suppressed by the smaller charge \Qhi — Qh { \- 
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Moreover, for the nucleus with large atomic number A > 20, it is usually dominated by spin- 
independent scattering [33]. Thus in the following discussion we only keep the contribution 



from spin-independent scattering. Note that the LDM-quark interaction is mediated by 
the dark gauge boson. Consequently, due to the conservation of the vector current, the 
sea quarks and the gluons will not contribute to the current operator f'fuf- As a result, 
the derivation of the effective theory is not only simplified, but also free of uncertainty 
from considerations like spin or strangeness content of the nucleon. This implies that the 
effective U(l)x charge of the nucleon, f pn , only receives contribution from its constituent 
quark |33|]. So we have 

f n = b u + 2b d , f p = 2b u + b d , (3.23) 

where b Uj d = dgy cos 2 QwQu,d with Q U: d being the QED-charge of u, d quarks. Thus at 
the leading order only the QED-charged proton carries a tiny U{\)x charge. The vector 
interaction only mediates the spin-independent interaction between dark matter and the 
nucleon and thus the LDM-nucleon scattering cross section can be added coherently to 
give the total LDM-nucleus cross section. That means that given the LDM-nucleon cross 
section a p , the LDM-nucleus cross section is proportional to (Zf p + (^4 — Z)f n ) 2 o p with 
Z and A being respectively the proton and atomic numbers of the nucleus. 

In practice, the four-fermion effective interaction is enough for the calculation of the 
LDM-nucleus cross section since, for each scattering by exchanging a dark gauge boson, the 
typical transferred momentum is (p± — P2) 2 = \q\ 2 = 2^ 2 v f 2 (l — cos#) ~ O(10 _4 )GeV 2 <C 
m 2 x ~ 4 x 10~ 2 GeV 2 . The reduced mass is /i/v = m^rriLi '(mi + m N) — tul for large 
nucleus like Ge with mass mx ~ 73 GeV. Thus the X boson can be integrated out. 

From Eq. (|3.22j ) we calculate the scalar spin-independent scattering cross section be- 



tween the LDM and the proton. It is given by 



N4 (Qngx) 2 



*p w T- 4 # ( 3 - 24 ) 
4 7T m % 



with N being the internal index from SU(N). This result is valid only in the non-relativistic 
limit (zero momentum transfer) . Note that in most previous studies the results are usually 
displayed on the plane of DM mass versus the DM-nucleon scattering cross section by 
setting f p = f n . But in our model, /„ w 0, and thus o p should be re-scaled as a p 
when compared with data, where w is the fraction of LDM in the total DM. Then we have 



4 9lQ 2 H h 9 2 x 

I — / ~~A 4 " ^ uav W 



/ Z 2 \ N m% glQ 2 Hh 9 2 x a2 



6 2 cos 9 



4 



72 



( z2 \ ... n„-,n-4o„ AT f0.2GeVy fQ H g x \ , u , 2 



, , , , x 1.2 x 10- 4U x N cm'. (3.25) 

\A 2 vu J \ m x J V °- 02 / V 2xl0_5 
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The cross section is independent of the LDM mass. Different experiments have different 
values for the ratio (Z/A) 2 . For a LDM with a mass of ~ 7 GeV, this cross section is just 
at the right order, according to the analysis in [1C]. 



4. Conclusion 

Both the cosmic ray anomalies observed by PAMELA and Fermi-LAT experiments and 
the possible events from direct detections like CoGeNT and CDMS II experiments may 
indicate the existence of dark matters. But the former points to a heavy dark matter 
at TeV sacle, while the later favors a light dark matter with a mass of several GeV. 
Meanwhile, the Sommerfeld enhancement may imply a dark U(l)x gauge boson with a sub- 
GeV mass. In light of these apparently different mass scales, we in this work constructed a 
simple and coherent framework with GMSB, based on the NMSSM extended with a light 
U(l)x sector and a heavy dark matter sector. By coupling the heavy dark matter directly 
to the U(l)x -charged messengers in the hidden sector, our framework has the following 
intriguing features: 

(1) The kinetic mixing 9 ~ 1CP 5 is obtained after integrating out the messengers with 
small doublet-triplet splitting. 

(2) A large negative mass-square m| for the NMSSM singlet S is generated at the mes- 
senger scale M. 

(3) The dark [/(l)xis spontaneously broken at the GeV-scale, while the dark gauge 
boson can have a sub-GeV mass by assuming a weakly U(l)x -charged Higgs field. 

(4) A GeV-scale light dark matter with the required interaction strength with quarks is 
provided. And its normal relic density can be generated by the ONLSP late decay. 
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A. Kinetic Mixing and Dark- Visible Interaction 

The messengers in the hidden sector are charged under U(1)y X U(l)x- After the SU(5)— 
gauge group is broken, these fields generate the kinetic mixing between the two groups. The 
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effective theory (for the bosonic gauge part only, with fermionic part obtained similarly) 
below the messenger threshold is given by 

C gauge = -\F$ v F Ylxv - \F%F X » V + °-F^F x , v . (A.l) 

Setting the initial value ft = at Mqut, then we get the value at the messenger scale M 
through one- loop RGE [34| 



rift 

16n2 % ^/Z^iQyQx, (A.2) 

where I runs over all the superfields with U(1)y x U(1)x charge. After electroweak sym- 
metry breaking, the three gauge bosons Z'„, B' and X'^ mix with each other through 
kinematic terms and mass terms. Using the convention in jnj], after eliminating the mix- 
ing and working in the mass eigenstate basis (Z^, A^, X^), the interactions between the 
gauge boson and the current at the leading order of 9 are described by 

Ccaupling D 9X„ (cos9 w J? m + 0(m 2 x /m 2 z )J%) + 9Z^ (- sw6 w J x + 0{m x /m 2 z )J^) 

+ 9 (bJ x + 0{m x /M x )XJ B ^ , (A.3) 

where the bosonic gauge current J^ mWZ an< ^ are defined as usual, while its fer- 
monic counterpart, i.e., the supercurrents, are defined as Jx = gx J2iQx,id\di, Jb = 
9Y Yli Qy,iv\vi with di/vi denoting any dark/visible fermions. The X^J^m accounts for the 
cosmic ray anomaly after HDM decays/annihilates to the dark sector. And BJx provides 
the U(1)y gaugino interaction with the dark states. 

B. Soft Terms from HDM-Messenger Direct Couplings 

If some fields feel the SUSY-breaking via direct couplings to the messengers, they will give 
new contribution to soft terms controlled by Yukawa couplings. The soft terms can be 
extracted by using wave- function renormalization method [22|. According to this method, 



we need to know the discontinuity of anomalous dimensions and the beta function across 
the messenger threshold [E4J. Above the messenger scale we have 



~<H h =-^- 2 {-2\l + AgxQ 2 H h ) 



lS h =T7~ 2 (-2X1 ~ 6A^ - 4X1 + ^\Q\ h ) 



16vr 2 
1 

lfr^ 2 

is ~ (-4k 2 - 2A 2 - 2A 2 ) , (B.l) 
2A 2 



_ 2X\ 
Xt 16vr 2 



5X1 + X 2 h + 2X1 ~ 2 9 \{Ql +Q 2 Hh ) ~ - Q 9l ~ yfll) , 



o\2 

Px D (3A 2 , + A 2 + 4X1 ~ 2g 2 x(Q 2 fl +Q 2 H h )-9l- ^l) ■ (B.2) 



-17- 



In the calculation, X is taken as a non-propagating background. Across the messenger 
threshold, the various discontinuity is given by 

A 7Hh =0, A 7Rh = -^(34 + 2A?,), A 7s = 0, 

9\2 

a /5a, = I ^(3A^ + 2Al ) ), A/3 A = A/3 K = 0, 

where c = (5/3, 1, 1) and n is the pair of (5, 5) messengers. Using this result, we calculate 
the following soft terms at the messenger boundary 

m l = - (1^)2 (3A| + 2X1) j^, (B.4) 
' (3A| + 2A 2 D )^, (B.5) 



16tt 2 v j M 

2 2,2 



k = (16^)2 t 8A ^ + 15A t + 12A|A|, - 16^ 



-6 32 2 A^ - 2g\ Qa| + A 2 , ) - 2 g |(Q 2 i + Q^)(3A| + 2A 2 ,) 



M 2 



>G- 



(B.6) 



(B.7) 



where m/Q = (Qn h /Q f 1 ) 2 m 2 ~/(2nx) (with dark gaugino mass m x = 2nx {&^r jj and 
nx = 5 in our paper) is the pure U(l)x GMSB contribution. And similar formula applies 
to the gauge mediation contribution to the soft mass term of other fields by replacing Qn h 
with corresponding charge. 

C. One-Loop RGEs of Some Soft Terms 

Here we present some important one-loop RGEs for all the soft terms in the dark sector, 
m|, A\ and A K in the NMSSM, and some Yukawa couplings and dark gauge couplings. In 
general, they take the form: 

dY 1 dA I dm 2 1 

where t = log(Q/Qo) with Qq being the boundary energy scale and Q the running scale. 
Following a general calculation in |35|], the RGEs for the soft terms in the dark sector are 
given by 

Pa Xi = 9A 2 A Ai +g 2 x Q 2 Hl {8m ^ - 4A Xl ) , (C.2) 
P m * H = 2NX 2 (m 2 Hi + 771% + m 2 H + A%) - 8g x Q 2 H m 2 ~ + 2Q H g 2 x S x , (C.3) 
m * Bt = 2A 2 (m 2 Hi + m\ x + m 2 H + A 2 ; ) - 8g 2 x Q 2 i m\ + 2Q Hl g 2 x S x , (C.4) 
P m 2 R = 2A 2 (m 2 Hi + m% + m 2 H + A%) - 8g 2 x Q%m\ + 2Q Sl g 2 x S x , (C.5) 
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where Sx is defined in the text. For the modified NMSSM soft parameters, the new RGEs are 
given by 

Pml = 4A 2 {m 2 Hu + m 2 Hd + m 2 + A 2 ) + An 2 (3m 2 + A 2 K ) 

+2X 2 (m 2 s + m 2 Hh + m\ h + A\ h ) , (C.6) 
p Ax = 8A 2 A A + 6h 2 A t + 6h 2 A b + 2h 2 T A T + 4n 2 A K + 2g 2 1 M 1 + 6g 2 M 2 + 2X 2 A Xh (C.7) 
Ak = 12k 2 A K + 12X 2 A X + 2X 2 A Xh , (C.8) 

Finally, for the new Yukawa couplings and the dark gauge coupling, their RGEs are given 
by 





= 2\\ (2k 2 + 3A^ + A 2 - A(Q Hh g x ) 2 ) (below M), 


(C.9) 




= 2A 2 (3A 2 - 4(Q Hh g x ) 2 ) (below M), 


(CIO) 




= 2g\ (lOQ 2 fl + 2Q 2 Hh + Q\ + N(Q\ + Q 2 Hi )) (above M), 


(C.ll) 




= fi,T (z£,t + ^T- yffl - ^9l) (above M), 


(C.12) 


%,£> 


= ii,D Uil, D + >?D- ~ ^gl) (above M). 


(C.13) 



At the scale Mqut, the unified value = £i ; t is assumed. 
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